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Abstract. Cellular automata have turned out to be important tools 
for the simulation of traffic flow. They are designed for an efficient im- 
pletmentation on the computer, but hard to treat analytically. Here we 
discuss several approaches for an analytical description of the Nagel- 
Schreckenberg (NaSch) model and its variants. These methods yield the 
exact solution for the special case v max = 1 of the NaSch model and are 
good approximations for higher values of the velocity (v ma x > 1). We 
discuss the validity of these approximations and the conclusions for the 
underlying physics that can be drawn from the success or failure of the 
approximation. 



1 Introduction 

Cellular automata (CA) do not only serve as simple model systems for the in- 
vestigation of problems in statistical mechanics, but they also have numerous 
applications for 'real' problems jjj. Therefore it is not suprising that in recent 
years CA have become quite popular for the simulation of traffic flow (see e.g. 

§)■ 

Here we do not want to discuss the question whether CA are suitable models 
for the description of traffic. We try to provide a tool box for the determination 
of the basic properties of the stationary state using analytical methods. These 
methods can also be applied to other CA models. 

CA are - by design - ideal for large-scale computer simulations. On the other 
hand, analytical approaches for the description of CA are notoriously difficult. 
This is mainly due to the discreteness and the use of a parallel updating scheme 
(which introduces 'non- locality' into the dynamics). In addition, these models 
arc defined through dynamical rules (e.g. transition probabilities) and usually 
one does not have a 'Hamiltonian' description. Therefore standard methods are 
not applicable. Furthermore, one has to deal with systems far from equilibrium 
which do not satisfy the detailed balance condition. 

However, there is a need for exact solutions or, at least, for good approxima- 
tions. These results as well as other exact statements may help to greatly reduce 
the need for computer resources. The interpretation of simulation data is often 
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difficult because of the 'numerical noise'. Even in the cases where an exact solu- 
tion is not possible, a combination of analytical and numerical methods might 
provide better insights. 

Here we will present analytical approaches which allow to solve the special 
case v m ax = 1 of the Nagel-Schreckenberg model || exactly. For v max > 1 these 
methods are only approximations, but they become exact in certain limits. We 
will discuss the quality of these approximations and how relevant information 
about the underlying physics can be obtained. 

2 Nagel-Schreckenberg model 

The Nagel-Schreckenberg (NaSch) model || is a probabilistic cellular automaton. 
Space and time are discrete and hence also the velocities. The road is divided into 
cells of length 7.5 m. Each cell can either be empty or occupied by just one car. 
The state of car j (j = 1, . . . , N) is characterised by an internal parameter vj 
(vj = 0, 1, . . . , v max ), the momentary velocity of the vehicle. In order to obtain 
the state of the system at time t + 1 from the state at time t, one has to apply 
the following four rules to all cars at the same time (parallel dynamics): 

Rl Acceleration: Vj — > v'a = Vj + 1 (only if Vj < v max ) 
R2 Braking: vj — > v'a = dj (for dj > Vj) 

R3 Randomization: v'a v'J = v'a — 1 with probability p (if > 0) 
R4 Driving: car j moves v'J cells. 

Here dj denotes the number of empty cells in front of car j, i.e. the gap or 
headway. One timestep t — > t + 1 corresponds to approximately 1 sec in real time 
H . For simplicity we will only consider periodic boundary conditions so that the 
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Fig. 1. Configuration in the Nagel-Schreckenberg model. The number in the 
upper right corner is the velocity of the car. 

number of cars is conserved^]. The maximum velocity v max can be interpreted 
as a speed limit and is therefore taken to be identically for all cars. 

1 A discussion of the effects of open boundary conditions can be found in the contri- 
bution by N. Rajewsky fl]. J. Krug fl and fj investigate the effects of disorder. 
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The four steps have simple interpretations. Step Rl means that every driver 
wants to drive as fast as possible or allowed. Step R2 avoids crashes between the 
vehicles. The randomization step R3 takes into account several effects, e.g. road 
conditions (e.g. slope, weather) or psychological effects (e.g. velocity fluctuations 
in free traffic). An important consequence of this step is the introduction of an 
asymmetry between acceleration and deceleration, e.g. overreactions at braking 
which are important for the occurance of phantom traffic jams. Finally, step R4 
is the actual motion of the vehicles. 

The NaSch model is a minimal model in the sense that all four steps R1-R4 
are necessary to reproduce the basic properties of real traffic. For more complex 
situations (e.g. 2-lane traffic) additional rules have to be formulated. 

In the following sections several methods for an analytical description of the 
NaSch model will be presented. 

3 Simple Mean-Field Theory 

The simplest analytical approach to the NaSch model is a (microscopic) mean- 
field (MF) theory. Here one considers the density c v (j, t) of cars with velocity v at 
site j and time t. In the MF approach, correlations between sites are completely 
neglected. 

For v m ax = 1 the MF equations for the stationary state (t — > oo) read jj: 

c Q = (c + pd)c, (1) 
ci = qcd (2) 

with d = 1 — c and q = I — p. The flow is simply given by /(c) = c\ = qc(l — c). 

For random-sequential dynamics the MF approach is known to be exact for 
Vmax — 1 || ■ For parallel dynamics, however, MF theory underestimates the 
flow considerably (see Fig. ||). This shows that correlations are important in this 
case. These correlations lead to an increase of the flow. Therefore one expects a 
particle-hole attraction (particle-particle repulsion), i.e. the probability to find 
an empty site in front of a car is enhanced compared to a completely random 
configuration. This picture will be confirmed by the exact solution. 

The MF equations for arbitrary v max > I are given in . One finds the same 
qualitative behaviour, i.e. MF theory underestimates the flow even more than 
for v max = 1. 

4 Garden of Eden States 

An important effect of the parallel dynamics is the existence of configurations 
which can not be reached dynamically § . These states are called Garden of Eden 
(GoE) states or paradisical states since they have no predecessor. An example 
for a GoE state is given in Fig. [| Note that the velocity is equal to the number 
of cells that the car moved in the previous timestep. For the configuration shown 
in Fig. this implies that the two cars must have occupied the same cell before 



Fig. 2. Fundamental diagram for v max — 1: Comparison of computer simulations 
(•) with the exact solution (continous line) and the mean-field result (broken 
line) . 

the last timestcpQ. Since this is forbidden in the NaSch model, the configuration 
shown can never be generated by the dynamics. 
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Fig. 3. A Garden of Eden state for the model with v max > 2. 



The simple mean-field theory presented in the previous section does not take 
into account the existence of GoE states. One can therefore hope that by elimi- 
nating all GoE states and applying mean-field theory in the reduced configura- 
tion space (paradisical mean-field, pMF) one will find a considerable improve- 
ment of the MF results. 

For v max = 1 all states containing the local configurations (0, 1) or (1, 1) are 
GoE states, i.e. the cell behind a car with velocity 1 must be empty. This only 
affects eq. (]l|) and the equations for pMF theory read: 

c = A/"(c +pd)c, (3) 
ci = Afqcd, (4) 

where the normalization J\f ensures Co + c\ = c and is given explicitly by Af = 
l 

Cg + d ' 

2 The reader should check whether Fig. hi depicts a GoE state or not ! 
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Fig. 4. Comparison of the fundamental diagrams obtain from MC simulations 
with the results of mean-field theory (MF) and mean-field theory without GoE 
states (paradisical mean-field, pMF). 

Solving (^) for c% by using cq = c — c±, one obtains c\ as a function of the 
density c. Since the flow is given by f(c) = ci, 

f(c,p) = l (l - Vl - 4g(l - c)c) (5) 

(with g = 1 — p) , we recover the exact solution for the case u maa; = 1 found 
first in || using the cluster approximation (see Section ||). Note that the flow 
depends on the density only via c(l — c), reflecting the particle-hole symmetry 
of the NaSch model for u mQ:E = 1 . 

For v max — 2 one has to take into account more GoE states ||. pMF is no 
longer exact, but it still leads to a considerable improvement of the MF results 
(see Fig. |). 

It is interesting to note that for random-sequential dynamics MF theory 
is exact whereas for parallel dynamics pMF is exact. This suggests that the 
important difference between random-sequential and parallel dynamics is the 
existence of Garden of Eden states in the latter. This is probably not only true 
for the NaSch model, but rather is a general property of CA models. 

5 Cluster Approximation 

The cluster approximation is a systematic improvement of MF theory which 
takes into account short-ranged correlations between the cells. In the n-cluster 
approximation a cluster of n neighbouring cells is treated exactly. The cluster 
is then coupled to the rest of the system in a self-consistent way. Related ap- 
proximations have already been used (under different names) for other models 
011111 • 

In order to simplify the description, we choose a slightly different update- 
ordering R2-R3-R4-R1 instead of R1-R2-R3-R4, i.e. we look at the system after 
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the acceleration step. Then there are no cars with v = and effectively we have 
to deal with one equation less. In the following we will use occupation variables 
rij where rij — 0, if cell j is empty, and rij — v, if cell j is occupied by a car with 
velocity v. 

If we denote the probability to find the system in a configuration (m , . . . ,til) 
by P(ni, . . . ,til) the 1-cluster approximation means a simple factorization 

L 

P(n 1 ,...,n L ) = l[P(n j ). (6) 

This is nothing but the mean-field theory of section || For the 2-cluster approx- 
imation one has a factorization of the form 

P{ni, . . . , til) oc P{n 1 ,n 2 )P(n 2 ,n 3 ) ■ ■ ■ P(n L -i, n L )P{n Ll rii). (7) 

The 3-cluster approximation is depicted graphically in Fig. ||. 



Fig. 5. The 3-cluster approximation. Shown are the central 3-cluster and all 
neighbouring 3-clusters with a nonvanishing overlap. 



In general, the master equation leads to (v max + 1)™ nonlinear equations in 
n-cluster approximation. Although this number can be reduced by using certain 
consistency conditions [0 , a solution is only feasible for relatively small cluster- 
sizes [Q. The quality of the approximation improves with increasing n and for 
n — * oo the n-cluster result becomes exact. 

However, for v max = 1 already the 2-cluster approximation is exact 
This is understandable from the results of the previous section since in this case 
the 2-cluster approximation is able to take into account all GoE states. The 
2-cluster probablities for the stationary state are given explicitly by 

P(0,0) = l-c-P(l,0), 

P(l,l) = c-P(l,0), (8) 
1 r 



P(1,0) = P(0, 1) = — [l - y/1 - 4 9 c(l - c) 

where again q = 1 — p. Using /(c) = qP(l,0), one obtains the exact result for 
the flow (see (||)). Note that the particle- hole attraction is apparent from these 
results, since P(0, 1) > P(0)P(1) = c(l - c). 
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density c 

Fig. 6. Comparison of simulation results with the n-cluster approximation (n = 
1, . . . ,5) for the fundamental diagram for v max = 2 and p = 1/2. 

For v max — 2 the fundamental diagrams obtained from the n-cluster ap- 
proximation (n — 1, ... ,5) are compared in Fig. || with results of Monte Carlo 
simulations. One can see a rapid convergence, already for n — 4 the difference 
between the simulation and the cluster result is extremely small. 

5.1 Correlation Length 

As an application of the cluster approximation we will compute density-density 
correlations for v max = 1 in the following. Using occupation numbers rij = 0, 1 
the density-density correlation function is defined by 

(nw r ) = nin r P(ni, ... , til) (9) 

where the prime indicates that the sum runs over all states with fixed particle 
number N = ri\ + ■ ■ ■ + Ul- 

In order to evaluate the sum in (^) , it is convenient to use a grand-canonical 
description. One introduces a fugazity z which controls the average number 
of particles (rij) and sums over all configurations in (^). Using the 2-cluster 
approximation (Q) - which is exact for v max = 1 — the correlation function is 
given by 

(nirir) = -^—'^2n 1 n r z N P(n 1 ,n2)P(n2,n 3 )---P(n L ^ 1 ,n L )P(n L ,n 1 ) (10) 

{ nj } 

with N = n j an d t ne normalization 

L 
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Introducing the transfer matrix 

- _ / P(0,0) V*P(0,1)\ , . 

F - ^P(1,0) zP(l,l) J 

this can be written succinctly as 

Z gc = Tr P L , (13) 



^- Tr (QP'- 1 QP L - r+1 ) . (14) 



^gc 

with the matrix Q(ni,ri2) = WiP(ni,ri2). 

The correlation length £ can be obtained from the asymptotic behaviour 
(r — > oo) of the correlation function 

(mn r ) - c 2 cx e~ r/ « (15) 

where c = (rij) is the (average) density of cars. £ is determined by the ratio of 
the eigenvalues A± of P (with |A+| > |A_|): 



£ _1 = In 



A_ 



Aj 



(16) 



The explicit expression for £ is rather lenghty. Therefore we just note that, for 
fixed p, £ is maximal for c = 1/2. £(c = 1/2) diverges only forp — * 0. In that case 
one finds that £(c = 1/2) cx p^ 1 ^ 2 . Monte Carlo simulations show that the same 
behaviour still occurs (at density c = - — — pj) for v max > 1 |lj,^5|. Therefore, 
the correlation function already gives an indication that the system is not critical 
for p > 0. The simulations show that there is no qualitative difference between 
the cases v max = 1 and v max > 1 as far as the phase transition is concerned |l5|| . 

The above results demonstrate that, although the 2-cluster approximation is 
exact, not all correlation functions are short-ranged or even of finite range. 



5.2 Further Applications 

We briefly mention other results which have been obtained analytically using the 
cluster approximation. In |l6| the distributions of gaps and the distance between 
jams have been calculated for v max = 1 using the 2-cluster approximation. For 
the gap distribution one recovers the (exact) result which will be derived in Sect. 
H (see eq. (]l8|)). For the calculation of the distribution of gaps between jams one 
defines every vehicle with velocity to be jammed. The distance between two 
jams is then given by the distance between a car with velocity and the next 
car with velocity 0. In |l7j the probability V(t) of a time headway t has been 
investigated. This quantity is defined in analogy to measurements on real traffic 
where a detector registers the time interval between the passing of consecutive 
cars. 
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6 Car-Oriented Mean-Field Theory 

The car-oriented mean-field (COMF) theory pl| is another possibility to take 
into account correlations in an analytical description^]. 

The central quantitiy in COMF is the probability P n (v) to find exactly n 
empty cells (i.e. a gap of size n) in front of a car with velocity v. In this way 
certain longer-ranged correlations are already taken into account. The essence 
of COMF is now to neglect correlations between the gaps. 

For v m ax — 1 the system of equations resulting from the master equation has 
the following form: 

P (t + l)=g(t) [P Q (t) + qP 1 (t)] , 

P x {t + 1) = g(t)P (t) + [qg(t)+pg(t)} P x (t) + qg(t)P 2 (t), (17) 
P n (t + 1) = pg{t)P n ^{t) + [qg(t)+pg(t)} P n {t) + qg(t)P n+1 (t), (n > 2) 

where g(t) = qJ2 n >i p n(t) = <?[! - Po{t)] (g(t) = 1 - g(t)) is the probability 
that a car moves (does not move) in the next timestep. As an example for 
the derivation of these equations we consider the equations for n > 2. Since 
the velocity difference of two cars is at most 1, a gap of n cells at time t + 1 
must have evolved form a gap of length n — 1, n, or n + 1 in the previous 
timestep. A headway of n — 1 cells evolves into a headway of n cells only if 
the first car moves (with probability g(t)) and the second car brakes in the 
randomization step (probability p), i.e. the total probability for this process is 
pg(t)P n -i(t). Similarly, the headway remains constant only if either both cars 
move (probability qg(t)) or both cars do not move (probability pg(t)). Finally, 
the headway is reduced by one, if only the second car moves (probability qg{t j). 

Although this infinite system of non- linear equations looks more difficult than 
those of the cluster approximation, a solution possible using generating functions 
[nsl. For v m ax — 1 one finds 



Po = 



1 
2q~c 



2qc-l + y/l - 4qc(l - c) 



P n = ^( j a - P t. ) (»>D, (18) 



p \P +p(l-P ) 

and for the flow f(c,p) = eg again the exact solution is reproduced. 

For v max = 2 one has two coupled systems of the type (|l8|), since one has to 
distinguish P n (v = 1) and P n (v — 2) fL8|| , This system can also be solved, but 
as for the other methods, one does not find an exact solution. 

In Fig. results for the distribution of headways are shown. For v max = 1 , 
the distribution has just one maximum at n = uq, with uq — 1 for small c 
and n = for large c. For higher velocities, however, there is a density regime 
where the headway distribution exhibits two local maxima Jl6| . The maximum 
at n = corresponds to jammed cars and the second maximum corresponds to 
free flowing cars. 

3 A similar method is discussed in Il9[. 
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Fig. 7. Distribution of headways for different densities for v max = 1 and p = 0.5 
(left) and v max = 5 and p = 0.5 (right). 



7 Cluster Approximation vs. COMF 



In this section the results of the cluster approximation and COMF in special lim- 
its are compared. As already mentioned, both methods yield the exact solution 
for Vmax = I- This is related to the fact that paradisical mean-field is exact in 
that case and both methods are able to identify all GoE states. For v max = 2 the 
situation is different. Here all three methods are only approximative. The quality 
of COMF is usually between those of the 2-cluster and 3-cluster approximations. 

In the limit p — > 0, COMF and the 3-cluster approximation become exact, in 
contrast to the 2-cluster approximation (see Fig. |||). Even for values of p w 0.1 
there is an excellent agreement between the fundumcntal diagrams obtained 
analytically and the Monte Carlo simulations. Since 'realistic' values of p are in 
the region p ~ 0.1 — 0.2 the approximations are indeed applicable in the relevant 
parameter regime. 

The limit p — > 1 is difficult to investigate numerically. Preliminary results 
show that COMF not exact in this limit. The cluster approximations give a 
much better agreement with Monte Carlo simulations, but it is not clear yet 
whether it becomes exact or not. 

In order to understand the differences between COMF and the cluster ap- 
proximation we calculate the distribution of jam sizes using these methods. Let 
C n be the probability to find a (compact) jam of length n, i.e. n consecutive 
occupied cells. C n is proportional to P(0|1)P(1|1) • • • P(1|1)P(1|1)P(1|0) in 2- 
cluster approximation. Here we have again used the different update-ordering 
R2-R3-R4-R1. P(ct|t) denotes the conditional probability P(a, r)/(J2 f P(er, f)). 
Note that the case n = 1 has to be treated separately, C\ oc J2l=T P(Q.\ v )P{n\Q)- 
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Cluster: v max = 2; p = 1/16 COMF: v max = 2; p = 0.1 




Fig. 8. Fundamental diagram for v max = 2: Comparison of simulation data (•) 
with the 2-cluster (• • • ) and 3-cluster results ( — ) for p = 1/16 (left) and with 
COMF (— ) for p = 0.1 (right). 

The number of jams is given by Afj = J2l=T an d one obtains 

= ^ f>(0>)P(<;iO) 

= ■^-P(0|l)P(l|ir- 2 P(l|l)P(l|0) (n > 2). (19) 

(2) 

Cn decays exponentially for n > 2, especially one has C n +i < C n . For the 
m-cluster approximation one can derive similar expressions. Now one finds an 
exponential decay for jam sizes n > m, i.e. small clusters may dominate (e.g. 
Cf } > Cf 5 is possible). 

In COMF, C n is proportional to (1 - P )P ■ • • P (l - P ) = (1 - P ) 2 P ™ _1 . 
The number of jams is propotional to 1 — P so that one finds 

C n = (l-P )P^- 1 . (20) 

This distribution is purely exponential and P„ > P n +i for all n. Therefore COMF 
is not able to describe clustering or phase separation, i.e. situations where jams 
with more than one car dominate. Clustering implies that there are correlations 
between gaps which are completely neglected in COMF. 

Since for p — > COMF is exact, there is no tendency towards clustering in 
this limit. On the hand, for p — * 1 COMF fails to give a good description of the 
phyisics since there is a strong tendency towards clustering or phase separation in 
that limit. Note that for the (deterministic!) case p = 1 the difference between 
Vmax = 1 and v max > 1 becomes most pronounced, since only for v max > 1 
metastable states exist [fl5i. 
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8 Other models 



In this section we discuss briefly two other models which have been investigated 
with the methods presented before. Both models are modified NaSch models. The 
first one is the T 2 model |2(],|2l| . It is a NaSch model with a modified acceleration 
rule. This so-called 'slow-to-start' rule is designed to mimic the delay of a car 
in restarting, i.e. due to a slow pick-up of the engine or loss of the driver's 
attention. We will discuss only the case v max — 1. Here the 'slow-to-start' rule 
has the important effect of breaking the 'particle-hole' symmetry of the NaSch 
model. 

Explicitly the 'slow-to-start' rule is given by: 

Rl' Acceleration with 'slow-to-start': A standing car will always accelerate to 
v = 1 if there are at least two empty sites in front; 

Is there only one empty site, it will accelerate only with probability 1 — p t . 

Fig. 9 shows the results for the fundamental diagram [^l| . Since the 'particle- 
hole' symmetry is broken, the flow is now maximal at a density c max < 1/2. 
Another interesting effect can be seen for larger values of pt- The flow vs. density 
relation has an inflection point and is no longer convex, i.e. the parameter p t 
controls the curvature of the fundamental diagram (Fig. 9). This can be seen 
most clearly for p t — 1. Here the flow f(c) vanishes for c > 1/2. This new phase 
transition to a completely jammed state exist only for p > 0. 




Fig. 9. Fundamental diagram for the T 2 model for different values of p and p t . 
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For c = 1/2 the stationary state consists of standing cars with exactly one 
empty site in between them (0.0.0.0. • ■ • where '0' denotes a standing car and 
'.' an empty site). No car can move due to the slow-to-start rule. For larger 
densities the stationary state is essentially of the same form, but larger clusters 
of standing cars will appear. In Monte Carlo simulations starting from a mega- 
jam the relaxation into this stationary state is extremely slow, even for 'large' 
values of p, e.g. p = 0.5. 

The COMF equations have two solutions Apart from the solution with 
vanishing flow for c > 1/2 there exists a second solution with /(c) > for 1/2 < 
c < 2/3 (dotted line in Fig. 9). Therefore one might expect the existence of a 
hysteresis effect in this model, which is indeed found in Monte Carlo simulations^. 
This is not surprising since the slow-to-start rule reduces the outflow from a 
megajam compared to the maximum flow of 'free traffic' which seems to be an 
important factor for the occurance of hysteresis. 

The second model has been introduced by Fukui and Ishibashi |2j| and can 
be considered as a NaSch model for "aggressive driving" . Again the rules are 
identical to the NaSch model, only the acceleration rule is changed: 

Rl" Acceleration: Every car accelerates to v max . 

For v max = 1 this does not change anything, but for higher velocities it leads to 
a considerable enhancement of the flow. 

Although the model is less realistic than the NaSch model, it is of interest 
due to its simplicity. An analytical description is much simpler than that of the 
NaSch model since after the acceleration step all cars have the same velocity 
Vmax- Therefore the analytical description is much simpler and might serve as 
testing ground for new methods p4fl. 




Fig. 10. Fundamental diagram for the Fukui-Ishibashi model: Comparison of 
simulation results (x) and COMF ( — ). 



4 A more detailed discussion can be found elsewhere in these proceedings |2: 
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9 Summary 

Although cellular automata are designed for efficient computer simulation stud- 
ies, an analytical description is possible, although difficult. We have present here 
four different methods which can be applied to CA models of traffic flow. The 
first approach, a simple mean-field theory for cell occupation numbers, is in- 
sufficient since the important correlations between neighbouring cells (e.g. the 
particle-hole attraction) are neglected. We therefore suggested three different 
improved mean-field theories. These approaches take into account certain cor- 
relations between the cells. The simplest method is the so-called 'paradisical 
mean-field' theory which is based on the observation that certain configurations 
(Garden of Eden states) can never be generated by the dynamical rules due 
to the use of parallel dynamics. The cluster approximation, on the other hand, 
treats clusters of a certain size exactly and couples them in a self-consistent way. 
Therefore short-ranged correlations are taken into account properly. In contrast, 
car-oriented mean-field theory is a true mean-field theory, but here one uses a 
different dynamical variable, namely the distance between consecutive cars. In 
that way, certain correlations between cells are taken into account. 

All three improved MF theories become exact for v max — 1 . For larger values 
of Vmax they are just approximations. In principle, the cluster approximation and 
COMF (in combination with a cluster approach) can be improved systematically. 
This is, however, very cumbersome. 

An interesting observation is that the qualitity of the approximation depends 
strongly on the value of p. This indicates that the physics changes with p, con- 
trary to common believe. Evidence for this scenario comes from the behaviour 
of the jam distribution and the phase transition. 

The methods presented here can also be used for other CA models. The 
investigation of the NaSch model has led to a better understanding of their 
advantages and limitations so that it is easier to choose the approach most 
suitable for a given problem. 

The question, how the stationary state is approached, is still an important 
open issue. Here, generalizations of the methods presented above can also be 
applied. This problem is currently under investigation |]25f and promises to yield 
new insights into the physics of the NaSch model. 

Acknowledgments I would like to thank Michael Schreckenberg for our fruitful 
collaboration which led to the development of the methods presented here. 
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